Competing many-body interactions in systems of trapped ions 
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We propose and theoretically analyse an experimental configuration in which lasers induce 3-spin 
interactions between trapped ions. By properly choosing the intensities and frequencies of the lasers, 
3-spin couplings may be dominant or comparable to 2-spin terms and magnetic fields. In this way, 
trapped ions can be used to study exotic quantum phases which do not have a counterpart in nature. 
We study the conditions for the validity of the effective 3-spin Hamiltonian, and predict qualitatively 
the quantum phase diagram of the system. 
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Most theoretical models in condensed matter physics 
rely on two-body interactions to describe a wide variety of 
phenomena, like for example, quantum phase transitions 
However, the presence of many-body interactions 
leads to exotic quantum effects, such as the existence 
of topologically ordered phases with anyonic excitations, 
which cannot be generally induced by pairwise couplings 
The beauty and complexity of those exotic models 
have motivated a good deal of recent theoretical work. 
Unfortunately, theory is still ahead of experimental im- 
plementations in this field, since many-body interactions 
are usually negligible in most systems found in nature. 

In this Letter, we show that 3-body effective spin in- 
teractions may be implemented with trapped ions, thus 
opening a new avenue of research beyond the traditional 
paradigm of condensed matter physics. The application 
of systems of trapped ions to the quantum simulation 0] 
of spin models has been theoretically analyzed, with a fo- 
cus on conventional models with 2-body interactions 0]. 
Furthermore, a recent proof-of-prinriple experiment has 
confirmed the validity of this idea [5]. By studying the 
implementation of 3-spin couplings, we show here that 
trapped ions can also be used to engineer new quantum 
states of matter, with a phenomenology that does not 
have a counterpart in nature. Recently, 3-body inter- 
actions have been theoretically analyzed in systems of 
ultracold atoms in triangular lattices [6] and cold polar 
molecules [7j. Our proposal would benefit from the ad- 
vantages of experiments with trapped ions, such as mea- 
surement and manipulation at the single particle level. 

To sketch our idea, let us consider a system of spins 
and a spin operator cr z = J2j a j • The spins are cou- 
pled to a bosonic mode (Ho = Sa^a) by both a linear, 
Hd = Fa z (a + at), and a quadratic squeezing term, 



H s = Mcr z (a 2 



J ). After the adiabatic elimination 



of the mode, the lowest order energy corrections in- 
clude 2-body (~ (F 2 /S)a z a z ,(M 2 /S)a z a z ), and 3-body 
(~ (MF 2 / 5 2 )cr z cr k erf ) interactions. In our scheme, this 
idea is realized in a system of ions in a linear array of 
microtraps [1, 0] . The internal states of the ions play the 
role of the effective spins. Two off-resonant laser beams 



are tuned such that they induce linear and squeezing 
couplings between internal states and the collective vi- 
brational modes of the chain. A third off-resonant beam 
induces an auxiliary linear coupling to partially cancel 2- 
body interactions, and thus, to tune the relative strength 
of the 2- and 3-spin couplings. In this way, experiments 
may reveal quantum effects beyond the usual pairwise- 
induced correlations, such as a quantum trier it ical point, 
and a quantum phase with a 3-body order parameter. 
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FIG. 1: Linear array of microtraps subjected to three laser 
beams along the axial and radial directions. 

We start by considering TV ions of mass m and charge 
e, confined along a one-dimensional array of microtraps, 
with lattice spacing a and trapping frequencies uo a with 
a = x,y,z (see fig. [1]). We assume that the ions have 
two internal hyperfine ground state levels (||), | j)), and 
their motion can be accurately described in terms of col- 
lective quantized vibrations (i.e. phonons). Therefore, 
the system Hamiltonian (h = 1) becomes 
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Here, h is an effective magnetic field induced by a laser or 
microwave field coupled to the internal transition, crj are 
the Pauli matrices corresponding to each ion, a^ na (a na ) 
are the phonon creation (annihilation) operators, and 
l] na = u a (l -\- c a f3 a V n ) 2 denotes the n-th normal mode 
frequency along the a-axis. These frequencies are ob- 
tained by diagonalizing V n = ^2 jk Mj n VjkMkn, where 
Vjk = ]jTfcp( 1 - 8jk) ~ Ez/j ]j=i\s S 3k is the Coulomb 
interaction in the limit of small vibrations, A4j n are the 
normal mode wavefunctions, and (3 a = e 2 /47reoracj^a 3 , 
c x ,y = 1, c z — —2. Note that in the "stiff" limit ^ Q < 1, 



2 



the trapping potential is much larger than the Coulomb 
repulsion, and thus vibrational frequencies are restricted 
to a narrow band of wid th ^c a /3 a co a around the trapping 
frequencies co a (see figs. |2(a)] and |2(b)| ). 

The ion chain will be subjected to three spin-dependent 
dipole forces Hi = Hl* + Hl x + ^i>f> sucn that the 
Hamiltonian becomes H = Hq + Hj. Here, the most 
general definition of a dipole force along direction a is 



tion Hamiltonian, Hx(t) = Hd(t) + H s (t), with 
H d (t) = £ F* n {t)a\ x o* + Fjn(t)al z o* + h.c. 



H s (t) 
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h.c, 



(3) 



Hra = 



2 



N 



where we have introduced the coupling strengths 
Ff n (t) = i^Lx^sincjLxt + Q L ^ 2 sin u; L ^t) M jne^ , 
and = i(Q L *rin 3 sin LJ L zt)Mj n e [cf>jz for the con- 



i(k L ar ja —uj L at) 



-\-hc^\cr z (2) tributions arising from the first-sideband, and analo 



which is induced by a pair of lasers in a Raman configu- 
ration, such that uo is the detuning between the Raman 
beams, kj> is the difference between the laser wavevec- 
tors [lOj], and is the two-photon Rabi frequency 11]. 
The position of the ions in ([2|) can be expressed in terms 
of normal modes, r ja = r® a + V -^ JTt f/?.t 



), 



where r*- = xoe x +jae z are the ion equilibrium positions. 
Let us stress that in the stiff regime, /3 a <C 1, it is possible 
to tune the dipole forces to a certain sideband of every 
normal mode of the ion chain, and thus couple each spin 
to the whole ensemble of vibrational phonons. 
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(a) Radial modes 



(b) Axial modes 



FIG. 2: Phonon frequencies Q na of (a) radial and (b) ax- 
ial modes in the stiff limit /3 a <C 1, where the corrections 
due to the Coulomb energy ^c a /3 a oo a V n have opposite signs 
for transverse c x , y = 1 or longitudinal c z = — 2 phonons. 
The coupling of spin-dependent dipole forces to vibrational 
modes is also shown: (a) In the radial stiff limit, we can si- 
multaneously blue-detune every mode to the first-sideband 
wl* ^ O nx , or to the second-sideband ujl x > 2Q nx . (b) In 
the axial stiff limit, every mode shall be red-detuned to the 
first-sideband ujl^ < Qnz- 

The main effect of Hl* and Hl% is to induce condi- 
tional linear and squeezing terms, which will produce 2- 
and 3-spin couplings, respectively. This is achieved by 
tuning the dipole forces to the first and second blue side- 
bands > Q nx , u L * = 2u L * > 2Q nx , see fig. |2(a)[ ). 
Besides, the role of Hl* is to generate 2-body couplings 
with an opposite sign as those induced by i^z*, and 
thus to partially cancel 2-spin interactions in favour of 3- 
body effects. As discussed below, the optimal screening 
is achieved by red-detuning this force to the first axial- 
sideband (ojli ^ ^nz, see fig. |2(b)[ )). In the limit of re- 
solved sidebands = \ ft na — uol™ \ <C H na , and weak 
couplings <c H na , we get a time-dependent interac- 



gously M/ nm (t) = \ (n L * cos UL*t)r)* 2 r)^ 2 M jn M jr ne 
for the second-sideband. Let us note that these cou- 
plings are switched on adiabatically, and depend upon 
the laser intensities through Ql™, the Lamb-Dicke pa- 
rameters rj^ K = k^c* I \j2mVt notl and the relative phases 
<j/j K = fcc^r^, assuming for simplicity 0^ = 0J 2 — <\> x - 

To get the explicit form of the 2- and 3-spin couplings, 
we first note that Ho + Hd(t) describes a set of forced 
harmonic oscillators. If we consider that dipole forces 
are adiabatically switched on, and t ^> l/(Q na — Ql k ), 
the evolution operator corresponding to Ho + Hd(t) is 

U d (t) = e-^e-^e 5 , with 
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where A = ±1 also considers non-resonant terms, and 

S = £ t>L°J " h - c -' $n = i f ; drF? n (r)e in »« T . (5) 
^ ■ , Jo 
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Let us emphasize here that the effective 2-spin interac- 
tions mediated by transverse (J|j > 0), and longitudinal 
phonons (JJ fe < 0), have opposite signs that interfere and 
lead to a partial screening of the two-body interactions. 
In this way, the effects of a 3-body coupling, which arise 
as a consequence of the additional non-linearities intro- 
duced by H s (t), can be experimentally accessed. Note 
that we are extending here the canonical transformation 
used in [3, Qjj], to the time-dependent case, something 
that may find a broader application in the context of 
trapped ion quantum simulation. Working in the inter- 
action picture with respect toUd-, the transformed Hamil- 

i^err, where 



t (3) 

tonian becomes U ] d H s Ud = # eff 
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describes the effective 3-spin interactions, and 



nmjk 



(7) 



-24Vfe)+h.c. 
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are the residual non-resonant spin-phonon couplings [14J . 

The conditions under which these error terms become 
negligible, which are thoroughly described below, lead to 
a novel Hamiltonian in quantum magnetism 



Et(3) z z z 
J jkl< J 3 a k a l 



(8) 



j>k 



j>k>l 



where 2- and 3-body interactions contribute. We note 
here that naive scalings ~ ry 2 , ~ n 4 (rj ~ 0.2 
in experiments) imply that two-spin couplings shall hin- 
der 3-spin effects. However, as shown in Eq. (|4]), the 
contribution to the 2-body terms from H L i (J? k > 0) is 
partially cancelled by that from H L 3 (Jj k < 0). This ef- 
fect becomes optimal if the parameters of the axial dipole 
force are carefully chosen. In particular, k^ should be 
an integer multiple of 27r/a, such that the ions sit on po- 
sitions with the same relative phase of the dipole force. 
To achieve ~ the laser must be red-detuned, 

with an intensity £l L 3 and detuning S z < fulfilling 

(fii s fci 5 )AK,(j z ,)9 z ) = (^|^| T )/ 2 (^,^,^), 

(nl z kl*)h(uj z ,5 z ,(3 z ) = (^|^l T )/ 3 (^,^,^)+ (9) 



where the functions fj(ou a ,5 a ,/3 a ) are listed in [13[. Fur- 
thermore, to obtain short-ranged spin interacting mod- 
els, the radial stiffness parameter is restricted to /3 a ~ 
0A5 a /cu a , where S a = uo a — uoL a is the detuning with 
respect to the bare trapping frequency. Under all these 
experimental constraints, at reach with current technol- 
ogy, we come to expressions for the many-body couplings 



Jjk = J2& jk , J)l\ = I h (AjkAki + A fej A jZ + AjiAik) , 

(10) 

where we have introduced the dipolar scaling function 
Ajk = — k\ 3 . The strength of these couplings is 



r(3) _ 1 



I* 



3f3£cj£F£M x cost 



(11) 



where x = AAM^, S z , Pz)/ Px\&x\h(u x , 5 X , f3 x ). We 
have introduced the bare linear and quadratic intensities 
by F x = ^Ql*Vx, M x = \ Q L xtjI, and the bare Lamb- 
Dicke parameter rj x = fc^/ '\/2mcj x . Hence, the access to 
the wavelengths, detunings, and laser intensities, leads to 
the controllability of 2 and 3-body couplings. 

In order to check the viability of this proposal, we 
should carefully deal with the residual spin-phonon cou- 
pling in Eq. Q, which contains different non-resonant 
terms that contribute to the error with 0(M 2 /c) 2 ) [THI ] 
and O / S x ) . Besides, the canonical transforma- 

tion in Eq. J5j leads to an additional error that scales as 
(9(F 2 /£ 2 ) [16[. Hence, a feasible quantum simulation of 
competing many-body interactions requires the parame- 
ters Mr, F x to be small in comparison to the laser detun- 
ing S x . At this stage, we can check the viability of this 



proposal assuming the following available experimental 
parameters, which fulfill all the constraints above. Con- 
sidering trapping frequencies oj x ~ 10MHz, laser detun- 
ings 5 X ~ 1.25MHz, stiffness parameters f3 x ~ 0.05, and 
setting the dipole force intensities to F x w M x ~ 0.18 x , 
we readily obtain an effective model of competing 2- and 
3-spin interactions with | Js\ ~ J 2 ~ 0.6kHz and an error 
on the order of E ~ 10 -2 . Furthermore, the modification 
of these parameters within the constraints detailed above, 
allows to experimentally access different regimes where 
I J3 1 / J 2 1 and thus observe the consequences of a com- 
petition between the many-body interactions in full glory. 
Let us note that one can also go beyond the usual dipolar 
regime in Eq. ([TP]) by considering j3 x u x /5 x > 0.5. In this 
case, exotic long range interactions between distant spins 
arise and offer a exceptional playground where the effects 
of the range of interactions can be studied. Additionally, 
the couplings can be raised to | J3 1 ~ J2 ~ 1 — 10kHz. In 
order to achieve this interesting regime, one may relax 
the trapping frequencies or design the microtrap in such 
a way that the ion equilibrium distance is lowered. 
The ability to independently tune the couplings 
^3? h) offers the opportunity to study novel quantum 
phases of interacting spins. To get a qualitative picture, 
we consider the Hamiltonian where only the nearest- 
neighbour terms of the dipole couplings are kept 



H 



J2 y^v 



+2 ■ 



■*£°?- (12) 

3 



The ground state is determined by the competition be- 
tween the different terms in ([T2|) : J2(> 0) induces anti- 
ferromagnetic (AF) order, Js(< 0) will be shown to in- 
duce a novel ferrimagnetic (F) phase, and h(> 0) encour- 
ages the system to lie in a disordered paramagnetic (P) 
regime, where spins are aligned along the x-direction. We 
start our analysis by considering the following two limits: 

(i) J3 = (2- spin quantum Ising model). This case 
is exactly solvable and shows a quantum phase transi- 
tion at the critical coupling J| = h, between the anti- 
ferromagnetic doubly degenerate ground state (J2 > h) 
\9af) £ {|U ••• t!),UT IT)}, and the paramagnetic 
phase \g P ) oc (g^d^j + (J2 < h) Q. To quantify 
the degree of anti- ferromagnetic order, we define the or- 
der parameter 0af(#) = j^T (dWj tf+i \d), which 
fulfills AF (g P ) = 0, and G AF (g AF ) = 1. 

(ii) J 2 = (3-spin quantum Ising model). In this case, 
3-spin interactions induce a novel quantum phase which 
can be fully characterisecharacterizedd by the order pa- 
rameter F (g) = 7^E^ri 2 (^l cr | cr |+i^I+2^)- The 3 " 
spin Ising model is no longer exactly solvable, but shows 
self-duality properties [17[, something that allows us to 
locate its critical point at the value |J|| = h. This 
point separates a phase with a four-fold ferrimagnetic 

state, \g F ) e {|TTT ••■ tTT>,IUI ••• IU),|m ••• ITI 

)» HIT ' • ' lit)} (J3 > h), from the paramagnetic phase 
(| Js\ < h). Remarkably, at Jf there is a phase transition 



which belongs to f the four-state Potts model universality 
class. Hence, 3-spin correlations induce an exotic critical 
behaviour different from the Ising universality class. 

From these limiting regimes, one gets a notion of the 
complexity of the model for general (J2, Js^h). In ad- 
dition to the critical points studied above, the system 
should also hold a quantum phase transition between 
the AF and F phases, as well as a trier it ical point, at 
which all the magnetic orders coexists (see fig. 3(a)[ ). This 
qualitative picture is supported by the finite-size numer- 
ical calculations presented in figs. |3(b)| and 3(c)| . Note 
that the order parameters, (9af, C?f> ma y be measured 
by detecting the photoluminescence from individual ions, 
something that amounts to a quantum measurement of 
<j z [3]. It is precisely this ability of performing highly 
accurate measurements at the single particle level, which 
allows us to characterize the full quantum phase diagram. 
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(b) AF order parameter (c) Ferri parameter 



FIG. 3: (a) Quantum phase diagram with P, AF and F phases 
that coexist in the tricritical point (in red). Order parameters 
as a function of the couplings J2, J3 for a chain with N — 15 
spins: (b) Anti-ferromagnetic order parameter (c) Ferrimag- 
netic order parameter. 

Let us briefly consider the implementation of these spin 
models using the more conventional linear Paul traps. Al- 
though the axial stiff limit cannot be achieved (i.e. the 
ion chain stability imposes f3 z > 1), it is still possible to 
devise effective two- and three-spin interactions. Focus- 
ing on the case of three ions, a similar procedure as that 
presented for microtraps, but considering a single mode 
in both radial directions a = x,y would yield 

H = J 2 (tTft7|+t7|t7j+t7j(7f) + J3^(7|(7|-h((7f +(7|+of ), 

where there is full access to the different couplings 
(^2,^3,^)- In this case, we could even switch off the 
two-body term J 2 = enhancing pure three-spin inter- 
actions and vice versa. In the same spirit as [5[ , one could 
perform a proof-of-principle experiment, where an initial 
separable paramagnetic state |P) = | adiabat- 
ically evolves towards an entangled state with different 
types of ordering. In case we tune J 2 < 0, J3 = 0, the evo- 
lution would generate GHZ states |GHZ) - | Ttt) + | III), 



whereas for J 3 < 0, J 2 > 0, one generates W- states 
\W) ~ I HI) + I HI) + I HI), studying thus the two 
non-equivalent classes of tripartite entanglement. Note 
that the Hamiltonians obtained with Paul traps are re- 
stricted to mean-field models, and microtraps should be 
used to implement strongly correlated systems. 

In conclusion, we have shown that a system of trapped 
ions can be used to explore the singular phenomenology 
of spin models with three-body interactions. We have 
made a realistic experimental proposal, at reach with 
current technology, to access the peculiar phase diagram 
of an effective Hamiltonian with three-body interactions, 
characterized by the appearance of exotic phases without 
a counterpart in usual condensed matter experiments. 
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